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ABSTRACT 

The galaxy intrinsic alignment is a severe challenge to precision cosmic shear measurement. We 
propose to s elf-calibrate the induc ed gravitational shear-galaxy intrinsic ellipticity correlation (the GI 
correlation. iHirata fc Seli"a3 12004) in weak lensing surveys with photometric redshift measurement. 
(1) We propose a method to extract the intrinsic ellipticity-galaxy density cross correlation (I-g) 
from the galaxy ellipticity-density measurement in the same redshift bin. (2) We also find a generic 
scaling relation to convert the extracted I-g correlation to the demanded GI correlation. We perform 
concept study under simplified conditions and demonstrate its capability to significantly reduce the 
GI contamination. We discuss the impact of various complexities on the two key ingredients of the 
self-calibration technique, namely the method to extract the I-g correlation and the scaling relation 
between the I-g and the GI correlation. We expect none of them is likely able to completely invalidate 
the proposed self-calibration technique. 

Subject headings: cosmology: gravitational lensing-theory: large scale structure 



1. INTRODUCTION 

of the most 
universe (IRefregieil 



Weak gravitational lensing is 

powe rfu l probes of t he da rk 

2001 lAlbrecht et all l200fil: iMunshi et al.l 120081 : 
Hoekstra fc JaiiJ l2008f l. It is rich in physics and 
contains tremendous information on dark matter, dark 
energy and the nature of gravity at large scales. Its 
modeling is relatively clean. At the multipole I < 2000, 
gravity is the dominant force shaping the weak lensing 
power spectrum while complicated gas physics only 
affect s the lensing power s pectrum at less than 1% 
level (iWhitel 120041: izhanin<:nox 200l iJing et al.ll2006l : 
iRudd et al.ll2008D . This makes the weak lensing precision 
modeling feasible, through high precision simulations. 
On the other hand, weak lensing has been measured 
with high significance. The most sophisticated method 
so far is to measure the cosmic shear, lensing induced 
galaxy shape di stortion. Af t er th e first detections in 
the year 2000 (iBacon et a,l.l 120001: iKaiser et al l 120001 : 
IVan Waerbeke et all 120001: IWittman et al.l \20m. data 
qua lity has been improved dramatically (e.g. IFu et al.l 

However, it is still challenging to perform precision 
lensing measurement. A potentially severe problem 
is the non-random galaxy intrinsic ellipticity (intrin- 
sic alignment). Its existence is supported by many 
evidences. Angular momentum of dark matter ha- 
los is influenced by the tidal force of the large scale 
structure. More importantly, mass accre tion is more 
prefe rentially along the filaments (e.g. IZhang et al.l 
l2009f ). This causes the dark matter halo elliptici- 
ties (and als o halo orientations) to be correlated over 



mentum vector, with t he halo shape , or wi t h each other 
(Hcvmans ct a ll 12004 [Yang et al.l 120(361: iKang et al.l 
i2007i : lOkumura et al.l 120081: IWang et al.ll2008D. t he re- 
sulting correlat ions in galaxy sp i n (iPen et al.l I200C ) 
and ellipticity (iBrown et al.l 120021: iMackev et al.l [ 2002 : 



[Hevmans et al."2004 'Hirata et al.'2004lLee fc PeniboOT : 
Okumura ct al. 2008; Schneider & Bridle 200^ are still 
measurable and contaminate cosmic shear measurement, 
especially for elliptical galaxies. 

The intrinsic alignment biases the cosmic shear mea- 
surement. It introduces an intrinsic ellipticity-intrinsic 
ellipticity correlation (the II correlation). Since this 
correlation only exists at small physical separation and 
hence small line-of-sight separation, it can be effec- 
tively eliminated b y the le nsing t omography technique 
dKing fc Schneided l2002l 120031 : iHevmans fc HeavensI 
2003), with moderate loss of cosmological information 
(|Takada fc Whitel[200l 

However, as pointed out by iHirata fc SeliakI (|2004) 
and then confirmed by ob s ervations ( fMandclbau m et all 
120061: IHirata et all l2007t lOkumura fc Jina .2009i) . the 
galaxy intrinsic ellipticity is also correlated with the 
large scale structure and thus induces the gravitational 
(cosmic) shear-intrinsic ellipticity correlation (the GI 
correlation). This contamination could bias the cos- 
mic shear measurement at 10% level and dark en- 
ergy constraints at more significant level ([Bridle fc Kingl 
[2003). For this reason, there are intensive efforts to 
correct the GI contamination (e.g. [Hirata & Sclia^ 
2004 iHcvmans ct al."2006': 'Joachimi fc Schneider 200§ 



20091: lOkumura fc Jina ,2009; Joachimi fc Bridle ,20091 : 
Shi et allboiol) ! 
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jCatelan et al.' 2001; .Crittenden et al.l 
Izhang et a l. 2009). Although galax- 
ies are not perfectly aligned with the halo angular mo- 
Electronic address: |pjzhang@shao.ac.cn| 



This paper proposed a new method to alleviate this 
problem. It is known that weak lensing surveys con- 
tains informatio n other than the ellipticity-e l hptici ty 
correlation (e.g. iHu fc JainI ([2004 ): iBernsteinI ([2009D ). 
These information not only helps improve cosmologi- 
cal constraints, but also helps reduce errors in lens- 
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ing measuremen t s, such as photometric r edshift errors 
(jSchneider et all l2006t IZhang et alJ I2OO90 . We show 
that these extra information allows for a promising self- 
calibration of the GI contamination. The self-calibration 
technique we propose relies on no assumptions on the 
intrinsic ellipticity nor expensive spectroscopic redshift 
measurements.^ It is thus applicable to ongoing or pro- 
posed lensing surveys like CFHTLS, DBS, Euclid, LSST 
and SNAP/JDEM. Under simplified conditions, we esti- 
mate the performance of this technique and its robust- 
ness against various sources of error. We find that it 
has the potential to reduce the GI contamination signif- 
icantly. 

The paper is organized as follows. In ^21 we explain 
the basic idea of the self-calibration technique. We in- 
vestigate the self-calibration performance in fJHl discuss 
extra sources of error in [2] and summarize in SJS] We 
leave some technical details in the appendices SjXl [B] & 

2. A SIMPLIFIED VERSION OF THE GI 
SELF-CALIBRATION 

In the section, we present the basic idea of the self- 
calibration technique. To highlight the key ingredients of 
this technique, we adopt a simplified picture and neglect 
many complexities until [JH 

2.1. The information budget in weak lensing surveys 

In lensing surve ys, we have several pieces of informa- 
tion (refer to iBernsteinl (|2009f) for a comprehensive re- 
view) . What relevant for the self-calibration technique is 
the shape, angular position and photometric redshift of 
each galaxy. Only galaxies sufficiently large and bright 
are suitable for cosmic shear measurement. To avoid 
possible sample bias, we will restrict the discussion to 
this sub-sample. We split galaxies into a set of photo- 
z bins according to their photo-z . The i-th bin has 
Zi — Azi/2 < z < Zi -\- Azi/2. Our convention is that, 
if « < J, then Zi < Zj. nf{z^) and ni{z) are the mean 
galaxy distribution over the i-th redshift bin, as a func- 
tion of photo-z z^ and true redshift z respectively. The 
two are related by the photo-z probability distribution 
function (PDF) p{z\z^). Intrinsic fluctuations in the 3D 
galaxy distribution, Sg, cause fluctuations in the galaxy 
surface density of a photo-z bin, (5e. This is one piece of 
information in weak lensing surveys. 

The galaxy shape, expressed in the term of ellipticity, 
measures the cosmic shear 71.2 induced by gravitational 
lensing. For each galaxy, this signal is overwhelmed by 
galaxy intrinsic ellipticity. If galaxy intrinsic ellipticity is 
randomly distributed, it only causes shot noise in the cos- 
mic shear measurement, which is straightforward to cor- 
rect. For this reason, we will not explicitly show this term 
in relevant equations, although we do take it into account 
in the error estimation. However, gravitational tidal 
force induces correlations in ellipticities of close galaxy 
pairs. For this, the measured shear 7^2 = 71.2 + -^i. 2, 

^ We still need suffici ently accurate plioto-z, which controls the 
I-g measurement error ^3.11 and the accuracy of the scaling rela- 
tion ( [133]). Since photo-z algorithm relies on calibration against 
spectroscopic samples (not necessarily of the same survey), in this 
sense, our self-calibration does rely on spectroscopic redshift mea- 
surements. 



where / denotes the correlated part of the galaxy in- 
trinsic ellipticity and 1 , 2 denote the two components of 
the cosmic shear and the intrinsic alignment. 71^2 de- 
scribe the shape distortion along the x-y axis and the 
direction of 45° away, respectively. 7 is equivalent to the 
lensing convergence a (in Fourier space, this relation is 
local). Thus we will work on k instead of 71.2. From 
the measured 7*, we obtain k" = k -\- 1. Here, / is the 
E m ode of Ii,2, ana. l ogous to k, which is th e E mode of 
71,2 (|Stebbinsl I199I: ICrittenden et al.ll2002[ ). Although 
cosmic shear, to a good approximation, does not have 
B-mode, the i ntrinsic alignment can have non- vanishing 
B-mode (e.g. iHirata fc Senaa[200l [Schneider fc Bridli 
[2OO9I) . This piece of information is useful to diagnose 
and hence calibrate the intrinsic alignment. However in 
the current paper, we will focus on the E-mode. 

K is the projection of the matter over-density along 
the line of sight. For a flat universe and under the Born 
approximation, the lensing convergence k of a galaxy 
(source) at redshift zq and direction 9 ()Refregieii I2003D 
is 

«(^) = 5nr{XL,e)WL{zL,ZG)dxL ■ (1) 

Here, 9 is the direction on the sky. (5m(xL, ^) is the mat- 
ter overdensity (lens) at direction 9 and comoving angu- 
lar distance xl = x(-Zl) to redshift z^. XG = xi^o) is 
the comoving angular diameter distance to the source. 
Both XL and xg ai'e in unit of c/Hq, where Hq is the 
present day Hubble constant. The lensing kernel 

WL{zL,ZG)^^a„,il + ZL)XL(l-—) . (2) 

2 V XgJ 

when zl < Zq and is zero otherwise. ^Im is the present 
day matter density in unit of the critical density. The 
lensing power spectrum is given by the following Limber 
equation, 

^^"^""W = 7 / ^™(^ = —^ZL)W,{zL)W,{zL)XLdXL ■ 
* Jo XL 

(3) 

Here, Wi is the lensing kern el W l averaged over the i-th 
redshift bin, defined by Eq. IA2I 

For two-point statistics, weak lensing surveys thus pro- 
vide three sets of correlations. Throughout this paper, 
we will work in the Fourier space (multipole i space) and 
focus on the corresponding angular power spectra. The 
first one is the angular cross correlation power spectrum 
between galaxy ellipticity (k^) in the i-th redshift bin 

and the one in the j-th redshift bin, cjf{i). 

Cl-\i) = C,'^«(^) + C^'^(^) + Cjf (^) + C!/ii) . (4) 

Here, C^J^ is the angular cross correlation power spec- 
trum between quantity a in the z-th redshift bin and 
quantity /? in the j-th redshift bin. a,l3 — G,I,g, where 
the superscript (or subscript in denoting the redshift 
and distance) G denotes gravitational lensing (k), / the 
galaxy intrinsic alignment (non-random intrinsic elliptic- 
ity) and g the galaxy number density distribution in the 
corresponding redshift bin (fe). 

The ellipticity-ellipticity pair relevant to the current 
paper has i < j. Since Cj/' <^ C-j^ as long as the catas- 
trophic error is reasonably small, we have 

cS^i) Cf^'^t) -f Clf[t) + 2ClHl) when i < j .(5) 
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For bin size Az > 0.2, as long as the catastrophic photo- 
z errors are reasonably small, Cl^ of non-adjacent bins 
(i < j — 1) is in general negligible with respect to the GI 
correlation (e.g. [Schneider & Bridle 2009), because the 
II correlation only exists at small line-of-sight separation. 
However, for adjacent bins {i = j — 1 ), the II correlation 
can be non-negli gible for Az - 0.2 (jSchneider fc Bridld 
l2009f ). The self-calibration technique proposed in the 
current paper is not able to co rrect for the II corre- 
lation, for which other methods (jJoachimi fc Schneideil 
[200a l2009t fZhandlMol) may be applicable. It only ap- 
plies to correct for the GI correlation C'^^ . We express 
the GI correlation as a fractional error to the lensing 
measurement, 

m^^§^y (6) 

flj is the fractional GI contamination to the lensing mea- 
surement, flj {i < j) and //j, {i < k) are not indepen- 
dent, since both describe the intrinsic alignment in the 
i-th redshift bin and thus 



(7) 



However, for its clear meaning as a fractional error in the 
lensing measurement, and for uncertainties in the intrin- 
sic alignment modeling, we adopt it, instead of C^j^ it- 
self, to express the GI contamination throughout the pa- 
per. The measured GI correlation is an anti-correlation 
ifij < 0)j because the lensing induced shear is tangential 
to the gradient of the gravitational potential while the 
intrinsic shear is parallel to the gradient. Throughout 
the paper, we often neglect this negative sign, since we 
will work under the limit f^j <C 1 and thus its sign does 
not affect our error analysis. Our self-calibration tech- 
nique works in principle for any value of /^^ . The results 
presented in this paper can be extended to other values 
of flj straightforwardly. 

The second correlation is between the galaxy density 
in the i-th redshift bin and the galaxy ellipticity 

(2) 

(k'*) in the j-th redshift bin, Clj . Galaxy-galaxy lensing 

in general focuses on C^-J^ (i < j). IZhang et al.l ()2009f ) 

showed that the measurement C^^ {i > j) contains valu- 
able information on photo-z outliers. The one relevant 
for our GI self-calibration is 



^gG 



gi I 



(8) 



The Cfl term contains here clearly contains valuable in- 
formation on the intrinsic alignment. The two terms on 
the right hand side has different dependences on photo-z 
error. Larger photo-z errors tend to increase C^*^ and 
decrease C^. 

The third set of cross correlation is between galaxy 
density ((5s) in the i-th redshift bin and j-th redshift 



bin, . It has been shown t hat Q,- {i ^ j ) is a 
sensitive probe of photo-z outliers ([Schneider et al.|[2006l : 
IZhang et aLll2009[ ). Our GI self-calibration requires the 
measurement 



(3) 



(9) 
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Fig. 1. — A schematic description of the GI self-calibration tech- 
nique. Here, G denotes gravitational lensing, / the non-random 
galaxy intrinsic alignment and g the galaxy number density distri- 
bution in the corresponding redshift bin. Here, galaxies in the j'-th 
redshift bin have higher photo-z than those in the i-th redshift 
bin. The GI contamination C^f^*^ is expressed as fractional error 

in lensing power spectrum , namely, f^j = /C^^ . The 

self-calibration operates when /^^ is bigger than certain threshold 
jthrcsh Depending on the fiducial value of ff^, residual GI con- 
tamination from different sources dominates, which we highlig ht as 
bold lines in lower part of the figure. Refer to Eg. [6l [T9l [211 fc [24l 
for definitions of corresponding variables. 



Our self-calibration aims to estimate and eliminate the 
GI contamination Cj{p in Eq. [5] from the measurements 

c'iP (Eq. [H]) and c|f^ (Eq. ©, both are available in 
the same survey. This method is thus dubbed the GI 
self-calibration. 

For the self-calibration to work, flj{(.) must be suf- 
ficiently large for the band power C^f{£) at the cor- 
responding £ bin to be detected through the measure- 
ment C'if \e). We denote the threshold as (f). For 
brevity, we often neglect the argument £ in fij. When 
fij > fij^'^°^^, we can apply the self-calibration to reduce 
the GI contamination. The residual GI contamination 
after the self-calibration is expressed in as residual frac- 
tional error on the lensing measurement, in which A/^ 
denotes statistical error and Sfij denotes systematical er- 
ror. Thus the self-calibration performance is quantified 
by f-j^'^°^^^, A fij and 6 fij. The smaller these quantities 
are, the better the performance. We will numerically 
evaluate these quantities later. 



2.2. The self- calibration 
The starting point of the GI self-calibration is a simple 



scaling relation between Cj{p and C^/ that we find 



h{t) 



elm 



(10) 
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Here, Xi = x(^i): ^i(^) is the galaxy bias in this red- 
shift bin at the corresponding angular scale i. Wij is the 
weighted lensing kernel defined by 



dzL I dzG[WL{zL,ZG)ni{zL)nj{zG)] 



(11) 

Notice that rii is normalized such that Jg ni{z)dz = 1. 
Ai is an effective width of the i-th redshift bin, defined 
by 

(12) 



Refer to the appendix fjA]for the derivation. 

This scaling relation arises from the fact that, both GI 
and I-g cross correlations are determined by the matter- 
intrinsic alignment in the i-th redshift bin. The prefac- 
tors in Eq. [10] are simply the relative weighting between 
the two. In this sense, this scaling relation is rather 
generic. 

The basic procedure of the self-calibration, as illus- 
trated in Fig. [H is as follows. 

• Extract C,-/ from the measurement C^^ in the i-th 

(2) 

photo-z bin. This exercise is non-trivial, since C^^ 
actually measures the sum of C// and C^^ , and 
C^^ is often non-negligible due to relatively large 
photo-z error. We find a method to simultaneously 
measure the two without resorting to spect-z infor- 
mation. The idea will be elaborated in and 
the measurement error in C^f will be calculated in 

• Measure the galaxy bias bi{£) from the measure- 
ment c^pie) (gaj). 

• Calculate Cjp from the above measurements and 
Eq. [To] and then subtract it from Eq. [4] 



2.3. 



Measuring Cj, 

An obstacle in measuring Clf comes from the con- 
tamination C^^ (Eq. [8]). For spectroscopic sample, this 
contamination is straightforward to eliminate. We just 
throw away pairs where the redshift of the galaxy to mea- 
sure the shape is lower than the one to measure the num- 
ber density ( Mandelbaum ct al. 2006; Hirata et al. 2007; 
lOkumura fc Jing|l2009[ ). This method is robust, however 
limited to the spec-z sample. 

For photo-z sample, the above technique does not 
work, since the photo-z error is large. For it, the true 
galaxy distribution, even for photo-z bin size Az — 0, 
has relatively large width >^ 0.1(1 + z), no matter how 
narrow the photo-z bin is. In practice, the photo-z bin 
size is often > 0.2, which further increases the effective 
width. Thus, a galaxy in this photometric redshift bin 
has large chance to lens another galaxy in the same red- 
shift bin, with a non-negligible lensing weight. For this 
reason, C^^ may not be negligible comparing to C//, 
even if the photo-z bin size is infinitesimal. We have 
numerically compared C^^ with C j,^ calculated based 
on the intrinsic alignment model of [Schneider fc Bridld 
(|2009f ) along their fiducial parameters, and found that 
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Fig. 2.— Qi{e) = C^i^\s{i)/C^i^{i) describes the suppression of 
the galaxy-galaxy correlation in the same redshift bin after throw- 
ing away those pairs with source redshift (photo-z) higher than lens 
redshift (photo-z). Only when Q deviates significantly from unity, 
can we extract the I-g correlation from the shape-density corre- 
lation measurement in the same redshift bin. Q is nearly scale 
independent, since it is the ratio of two power spectra of similar 
shape. Overall, Q ~ 1/2 . It increases with redshift, as expected 
from larger photo-z rms error at higher redshift and hence larger 
effective redshift width. This figure is a key result to demonstrate 
the feasibility of the self-calibration technique. 



the two terms can indeed be comparable. For example, 
at a typical lensing angular scale £ = 10'^ and a typical 
redshift z = 1, C^iyclf is 30% for vanishing bin size 
and increases with bin size. Only when the redshift is 
sufficiently low or the redshift error and the bin size are 
both sufficiently small, may the G-g correlation be safely 
neglected. 

Nevertheless, the photo-z measurement contains use- 
ful information and allows us to measure C^f . Our 
method to separate C^^ from C^f relies on their dis- 
tinctive and predictable dependences on the relative po- 
sition of galaxy pair. Let's denote the redshift of one 
galaxy in the pair for the shape measurement as Zq and 
the other one for the number density measurement as 
Zg . The I-g correlation does not depend on the ordering 
along the line-of-sight, as long as the physical separa- 
tion is fixed. In another word, the I-g correlation for 
the pair with Zq > z^ is statistically identical to the 

z^\ fixes. On the 



pair with Zq < Zg, as long as \zq 

other hand, the G-g correlation cares about the ordering 
along the line-of-sight. The G-g correlation for the pair 
with Zq > z^ is statistically larger than the pair with 

Zq < Zg , due to the lensing geometry dependence. 

We can then construct two observables on the 
ellipticity-density correlation. (1) One is C^f\£), weight- 

ing all pairs equally. (2) The other is C^^ \s{^), in which 
we only count the cross correlation between those pairs 
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with Zq < Zg . The subscript "S"' denotes the corre- 
sponding quantities under this weighting. From the ar- 
gument above, we have C//(£) = C//|s(£). On the other 
hand, C^^ls(^) < C^.^ii), since those zg > z^ pairs 
that we disregard contribute more to the lensing-galaxy 
correlation. We denote the suppression by 



(13) 



Q is sensitive to photo-z errors. In general catastrophic 
errors drive Q towards 1. Q = 1 if the photo-z is com- 
pletely wrong and has no correlation with the true red- 
shift and Q = if the photo-z is 100% accurate. Usually 
< Q < 1. Q can be calculated from the galaxy rcdshift 
distribution (the appendix fjB]) and is thus in principle 
an observable too. From the two observables 



Gg, 



We obtain the solution to C^f , 



1 — r<^9 1 



Gg\ 



1 - Qie) 



(14) 



(15) 



For it to be non-singular, Q must deviate from unity 
(Q < 1). We numerically evaluate this quantity and find 
that, for the survey specifications presented in this paper, 
Q ^ 1/2 (Fig. [2]). The significant deviation of Q from 
unity is not a coincidence. It is in fact a manifestation 
that the lensing efficiency changes significantly across the 
redshift interval crp. We thus expect the C-f estimator 
(Eq. [15]) to be applicable in general. 

3. ERROR ESTIMATION 

Unless explicitly specified, we will target at LSST 
throughout this paper to estimate the performance of the 
self-calibration technique proposed above. LSST plans to 
cover half the sky (/sky = 0.5) and reach the survey depth 
of ^ 40 galaxies per arcmin^ . We adopt the galaxy num- 
ber density distribution as n{z) oc exp(— z/0.5), the 
rms shape error 7rms = 0.18-1- 0.042z and photo-z scat- 
ter dp — 0.05(1 + z). We split galaxies into photometric 
redshift bins with Azi = 0.2 centered at Zi = 0.2{i + 1) 

(i = 1,2, •••)•' 

3.1. Measurement error in C^f 

('2) (2) 

Both and | s have cosmic variance and shot 
noise errors, which propagate into C^f extracted from 
the estimator Eq. [151 The error estimation on C// is 
non-trivial since errors in C^^ and ^ | s are neither 

^ The choice of redshift bin is somewhat arbitrary. For example, 
we can include a lower redshift bin with Zi = 0.2 (0.1 < < 0.3). 
The self-calibration certainly works for this bin, since the I-g cross 
correlation in this bin is easier to extract due to weaker lensing 
signal and hence weaker G-g in this bin. The major reason that we 
do not include this bin is that the lensing signal in this bin is weak. 
Weak lensing signal may cause confusion on the performance of 
the self-calibration, due to our choice to express the GI correlation 
before and after the self-calibration as ratios with respect to the 
lensing signal. For example. Fig. |3] shows that Afij increases 
toward low redshifts. But it is an manifestation of weak lensing 
signal instead of poor performance of the self-calibration. 
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Fig. 3. — The threshold of the applicability of the self-calibration 
technique /'j"^'^'^'^ and the residual statistical uncertainty A/^^"^ 
(solid lines with data points). Refer to Eq. llQI for the definition of 
A/j'°' . The self-calibration applies for sufficiently large GI contam- 
ination (/^^ > /*'^°'''^) and reduces the fractional error from f^^ to 
A/?j\ Here, the superscript (a) denotes the error from the 



measurement. Notice that ff^''^'^^ 



and both are insensi- 



tive to the fiducial /^^ . The cosmic variance in the lensing field and 
the random shape fluctuation sets the minimum fractional statis- 
tical error e^^'" in C^^'^ measurement, which we plot as dash lines. 
The numerical estimation shown is for redshift bins Zi = 0.2{i + 1) 
(i = 1, 2, ■ ■ ■) and multipole bin size A£ = 0.2^ in LSST. In general 

Af.'"' < e™'" and thus the residual error after the self-calibration 



is negligible. Since both A/^ 



(a) 



and e"J'" 



scale in similar ways, this 



conclusion holds in general for other lensing surveys. 



completely uncorrelated, nor completely correlated. The 
derivation is lengthy, so we leave it in the appendix iJCl 
The final result of the rms error AC^f in a bin with width 
Ai is 



(AC^) 



ig\2 _ 



1 



f<ggf< 

,r<99'^r<GG.N 
-t-O,-,- L/,-,- 



GG 



1 



1 



3(1 



1 



1 



) 
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Here, the superscript "N" denotes measurement noise 
such as shot noise in galaxy distribution and ran- 
dom shape noise. Cf/'^ — Airfsky/Ni and Cf^^'^ — 

47r/sky7rms/-^i, where Ni is the total number of galaxies 
in the i-th redshift bin. 

From Eq. [TH AC// is insensitive to the intrinsic align- 
ment, as long as the II correlation is sub-dominant with 
respect to the GG correlation. We will work at this limit. 
If the intrinsic alignment is too small, the measurement 
error AC// will be larger than . AC^f = C// hence 
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sets a threshold //^ . Combining Eq. [TOl and the defini- find that this is indeed the case in general Fig. [3] is the 
tion //■ = Clf jCf^ (Eq. [H), we obtain this threshold 



as 



/•thresh 
/ ij 



r<GG 



(17) 



jthrcsh j^g^g important meanings. (1) It describes 
the minimum intrinsic alignment (/^^ — that can 

be detected through our method with S/N=l. Thus 
it also defines the lower limit beyond which our self- 
calibration technique is no longer applicable. (2) It de- 
scribes the self-calibration accuracy resulting from C// 
measurement error. The measurement error AC// prop- 
agates into an error in C/^^ determination through Eq. 
[TUl and hence leaves a residual statistical error in the 
GG measurement. We denote this error as A/^^j*''. Since 

^fu ^/fh ^ ^Cu/Cl!^ combining Eq. M we have 



(18) 



Also, we find an important relation between the two, 




A/, 



(a) 



thresh 



(19) 



This relation can be understood in this way. When 
the intrinsic alignment fl^ > it can be inferred 

through measurement in C^f and hence be corrected. 

Since measurement C// has statistical error AC// (Eq. 
I16p . this correction is imperfect. The residual error in 
flj is set by the same AC// determining Jthrcsh^ 
have the relation Eq. [T21 If only our method is applied, 
we are only able to detect the intrinsic alignment with 
an amplitude f,lj > /'hrcsh ^^^^ correct it to a level of 

^yr(ti) ythresh 

Once the intrinsic alignment is sufficiently large to be 
detected (//,- > fij"''°^^), our self-calibration technique 
can detect and thus eliminate the GI contamination. It 
renders a systematical error in lensing measurement with 
amplitude f^^ into a statistical error with rms A/^^j^^ = 
jthresh ^ j/^ j^oi-i^g ^Yiat Af^f is insensitive to //^.. 

Numerical results on Afl""^ = /*'"°^^ are evaluated 
through Eq. [18] along with Eq. [16] and are shown in 
Fig. [3] For cornparison, we also plot the minimum G-G 
error va the measurement. It is the rms fluctuation 
induced by the cosmic variance in the G-G correlation 
and shot noise due to random galaxy ellipticities, in ide- 
alized case of no other error sources such as the intrinsic 
alignment. This minimum G-G error sets the ultimate 
lower limit of the fractional measurement error on C'^'^ 



a 



GG,2 



GGa 



(20) 

When Afj^j ^ is smaller than e™'", the residual error after 



^min 

^ij 



the self-calibration will then be negligible, meaning a self- 
calibration with little cosmological information loss. We 



forecast for LSST.^ Since both Af^f and e™'" scale in a 
similar way with respect to survey specifications such as 
the sky coverage and galaxy number density, this conclu- 
sion also holds for other lensing surveys such as CFHTLS, 
DES, Pan-Starrs, Euchd and JDEM. 

How can the GI contamination be corrected to an accu- 
racy even below the statistical limit of GG measurement? 
Equivalently, why is flj as small as the one shown in Fig. 
[3] detectable? This surprising result requires some expla- 
nation. The reason is that C// is amplified with respect 
to C/jG by a large factor l/WijA, - 0(10^) (Eq. [TO]). 
Thus a small GI contamination {f^j ^ 1) can still cause a 
detectable C//. This explains the small f^j''''^'^ = Af^f . 

3.2. Measuring the galaxy bias 

The second uncertainty in the self-calibration comes 
from the measurement of the galaxy bias, Abi{£). From 
Eq. 1101 this uncertainty induces a residual statistical 
error 

A/,^) = 4(A6,(^)/6.(^)) . (21) 

There are several possible ways to obtain hi such as 
combining galaxy-g alaxy lensing and ga laxy clustering 
measurements (e.g. iSheldon et al.l [20041) . c ombining 2 
point and 3-point galaxy clustering (e.g. (iGuo fc Jing 



2009)), fitting against the halo occupation (IZheng et al. 
2005|i and conditional luminosity function (jYang et al 
20031). or a nalvzing the counts-in-cells measurements 
( Szapudi fc Fan! 1200 41 . Alternatively, hi {£) can be in- 



ferred from the galaxy density-galaxy density correla- 
tion measurement alone, C'>f\£) = Cff{e) ~ bjC^'^ii). 
Here, C™™ is the matter angular power spectrum, with 
the same weighting as galaxies. Both uncertainties in the 
theoretical prediction of C™™ and measurement error in 

(3) 

affects the bi{l) measurement. Given a cosmology, 
one can evolve the matter power spectrum tightly con- 
strained at the recombination epoch by CMB experiment 
to low redshift and thus predict C™™. As long as the as- 
sociated uncertainty is smaller than 10%, the induced 
error will be sub-dominant to the systematical error dis- 
cussed later in i|3.3l This is an important issue for further 
investigation. 

On the other hand, the statistical error in bi {£) induced 
by the galaxy clustering measurement error is 



Ab,{e) 




1 



tAtU^y 



(22) 



This rough estimation suffices for the purpose of this 
paper, for which the reason will be come clear latter 
soon. This error is negligible for a number of reasons. 
(1) It is in general much smaller than the systematical 



error in the scaling relation, 5f^^ 



As will be shown 



0.1//-. On other other hand, 6^(^) can 



in general be measured with better than 10% accuracy. 

We do notice that in some cases especially when one of the 
photo bin is at low redshift, A/^^j"'' > e™"' (Fig- O, leading to 
non-negligible loss in cosmological information for relevant redshift 
bins. 
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For example, for LSST at ^ = 100 with M = 0.2^, 
Abi{£)/bi{£) ~ 1.6%. can be measured with higher 

accuracy toward smaller scales, until shot noise domi- 
nates. Thus Af^!^^ < Sf^f at relevant scales. (2) It 
is smaller than the minimum statistical error in 
measurement (Fig. [3]). First of all, galaxy cluster- 
ing measurement is in general more accurate than lens- 
ing measurement. Second, the impact of uncertainty in 
bi{£) on the self-calibration is modulated by a factor f-j 

{Aflf = fi^{Ah{£)/hm- Unless > 1, Aflf is 
suppresed. Thus we expect that the hi{t} induced er- 
ror Aff^^ is negligble even to the minimum statistical 
error in C^'^ measurement. From the above general ar- 
gument, this conclusion should hold for most, if not all, 
lensing surveys. We show one example of LSST. Even for 
a rather large fl^ = 1, uncertainty in bi{£) only causes a 

statistical error of Afj;^^ = 1.6% at £ = 100, negligible 
comparing to statistical uncertainties in Cf^'-^ measure- 
ment (Fig. [3]). The above conclusions are safe even if 
Eq. [21] underestimates the error in hi {£) by a factor of 
a few. This is the reason that we do not seek a more 
robust estimation on the measurement error in bi{£). 

From the above argument, the errors AJ^"^^ and 

^^if' from different sources and thus are inde- 

pendent to each other. The combined error is then 

(A/j^°^)2 -f {Afl'^^Y. Since the galaxy bias induced er- 
ror is likely sub-dominant to either the error source (a) 
or to the one will be discussed in ij3.3l we will neglect it 
for the rest of the paper. 

3.3. The accuracy of the Cfj^-C^f relation 

A key ingredient in the self-calibration is Eq. [TOl which 
links the observable C// to the GI contamination. How- 
ever, Eq. [10] is not exact and we quantify its accuracy 

by 



hmt^{£) 
W.,A,Cl!{£) 



1 



(23) 



also quantifies the induced residual systematic error, 

(24) 

We present a rough estimation by adopting a toy model 



AlAk, z) cx A2 (fc, z) cx A,2„(fc, z)(l + zf (25) 



with P = -1,0,1. Here, A^^, A^^ and A^ are 
the 3D matter-intrinsic alignment, galaxy-intrinsic align- 
ment cross correlation power spectrum (variance) and 
matter power spectrum, respectively. The accuracy of 
Eq. [TU] is not only affected by the scale dependence of 
corresponding 3D power spectra, but also their redshift 
evolution (the appendix \K\. Theoretical models of the 
intrinsic alignment (e.g. [Schneider fc Bridlil2009D show 
that the redshift evloution may not follow the evolution 
in the density field. For this reason, we add an extra 
redshift dependence (1 -|- zY in Eq. [5S] This recipe is 
completely arbitrary, but it helps to demonstrate the ro- 
bustness of Eq. [To] 
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Fig. 4. — €ij quantifies the accuracy of Eq. 1101 which links 
the observable C\? to the GI contamination C^fp. It thus quanti- 
fies a dominant systematical error of the self-calibration technique, 



5/> 



^iiSly Solid, dotted and dash lines represent three toy 



models of the intrinsic ellipticity evolution. Usually, Eq. llOl is accu- 
rate within 10%. However, for those adjacent bins with i, j = t + 1, 
due to stronger redshift dependence of the lensing kernel, Eq. 1101 
is least accurate and e^j is largest (right panel). 



Numerical results are shown in Fig. [3] We see that for 
most ij pairs, \ is less than 10%, meaning that we are 
able to suppress the GI contamination by a factor of 10 
or larger, if other errors are negligible. 

We notice that the largest inaccuracy of Eq. [TU] occurs 
for those adjacent bins (z, j = i -I- 1), which is 10%-20%. 
This is caused by the lensing geometry dependence. Eq. 
[TUI would be quite accurate if the integrand in Eq. lAll & 
IA3I varies slowly with redshift. However, since the i-th 
bin and the J = * + 1-th bin are close in redshift, the 
lensing kernel Wj{z) varies quickly over the zth redshift 
bin. This fast redshift evolution degrades the accuracy 
of the scaling relation and thus causes a larger \eij\. On 
the other hand, Wj{z) changes more slowly over other 
bins with z 7^ j — 1, since now the source-lens separation 
is larger. For this reason, Eq. [10] is more accurate for 
these bins. The self-calibration technique does not work 
excellently for adjacent bins, but a factor of 5 reduction 
in the GI contamination is still achievable. 

The scaling relation accuracy is senstive to the photo-z 
accuracy, eij cx tr^ ~ o-p -|- (Az)^/12. Here, a is the rms 
redshift dispersion in the corresponding redshift bin. One 
can obtain the above relation by perturbing Eq. lAll & 
I A3I around the median redshift and keeping up to the sec- 
ond order. Thus, if photo-z accuracy can be significantly 
improved, the accuracy of Eq. [TU] can be significantly 
improved too. For example, if <7p — 0.03(1 -I- z) instead 
of the fiducial value 0.05(1 -I- z), the scaling relation ac- 
curacy can be improved by a factor of ^ 2-3. This would 
allow us to suppress the GI contamination by a factor of 
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10-20 or even higher. The shape of the photo-z PDF also 
matters. Depending on which direction it is skewed, the 
scahng relation accuracy may be improved or degraded. 
In i )4.21 we will discuss the impact of catastrophic error, 
which presents as significant deviation from the adopted 
Gaussian PDF. 

The accuracy of the scaling relation may also be im- 
proved by better modeling, has (much) larger chance 
to be positive (Fig. 2]). This behavior is likely general, 
not limited to the intrinsic ellipticity toy models we in- 
vestigate. In deriving the scaling relation, both C^f and 
Clj^ are evaluated at the middle redshift Zi. The weight- 
ing function in C^f roughly peaks at Zi while the one 
in peaks at lower redshift, due to the monotonous 
decreasing of Wj{z) with redshift (until it vanishes). It 
is this imbalance causing the general behavior > 0. 
It also shed ligh on improving Clj^-C^f relation and re- 
duceing the associated systematical error: an interesting 
project for future investigation. 

3.4. General behavior of the residual errors 
Depending on the nature of the GI correlation, either 
the error in C^f measurement or the error in the scal- 
ing relation can dominate the error budget of the self- 
calibration, while the one from the galaxy bias is likely 
sub-dominant. Fig. [T] demonstrates the relative behav- 
ior of the three error sources in the self-calibration. The 
bold lines highlight the dominant error, as a function of 
the intrinsic alignment amplitude f^^. There are several 
regimes. 

• fij < fij^'^'^^^- The intrinsic alignment is too 
weak to be detected in the galaxy-lensing corre- 
lation. For this reason, the self-calibration tech- 
nique is not applicable. However, this usually 
also mean the intrinsic alignment is negligible in 
lensing-lensing measurement, comparing to the as- 
sociated minimum statistical error (Fig. [3]). Thus 
there is no need to correct for the GI contami- 
nation in this case. However, there are impor- 
tant exceptions to the above conclusion. For ex- 
ample, from Fig. |3l when one photo-z bin is at 
sufficiently low redshift, the GI contamination is 
undetectable by our method, but the systemat- 
ical error it induces is non-negligible. Further- 
more, our method is insensitive to the intrinsic 
alignment which is weakly correlated to the den- 
sity field. Such intrinsic alignment can cause large 
contamination to the lensing power spectrum Cff'-^ , 
but leaves no detectable feature in the ellipticity- 

(2) 

density measur ement C;.- . In these cases , other 
methods (e g. iJoachimi fc Schneided I2008L 120091 : 
IZhang|[2bl0( ) shall be applied to correct for the in- 
trinsic alignment. 



thresh 



The self-calibration begins to work. 



(1) flj < A/*j^'''="h/eij. The statistical error in- 
duced by Tg measurement uncertainty dominates. 



However, this residual error. A/, 



(a) 



A/, 



thresh 



usually negligible, comparing to the minimum sta- 



e™™. Fig. [3]). In this domain, the self-calibration 
technique is promising to work down to the statis- 
tical limit of lensing surveys. 

(2) fIj > A/*^j^^'^"i>/e,j . The systematical error aris- 
ing from the imperfect scaling relation domiantes. 
The fractional residual error in lensing-lensing mea- 
surement is Sflj^ = eijflj ~ 0.1/j^-. This error will 
be still sub-dominant to the lensing statistical fluc- 
tuation, if f,fj < efij™/eij. If not the case, the self- 
calibration can work to suppress the GI contami- 
nation by a factor of 10. Other complementary 
te chniques such as the nulling techniq ue proposed 
bv IJoachimi fc SchneideH (|2008l[2009l ) shah be ap- 
plied to further reduce the residual GI contamina- 
tion down to its statistical limit. 

4. OTHER SOURCES OF ERROR 

There are other sources of error, beyond the three ones 
discussed above. We discuss qualitatively on magnifi- 
cation bias, catastropic photo-z error, stochastic galaxy 
bias and cosmological uncertainties. Based on simplified 
estimations, we conclude that none of them can com- 
pletely invalidate the self-calibration technique. Quanti- 
tative and comprehensive evaluation of all these errors, 
including the ones in previous section, will be postpone 
elsewhere. 

4.1. Magnification bias 

Gravitational lensing not only distorts the shape of 
galaxies, but also alter their spatial distribution and in- 
duces magnific ation bias, or cosmic magnification (e.g. 
IScranton et al.l ()2005i) and references therein) . It changes 
the observed galaxy overdensity to 5g = 5g + g{F)n. The 
function g{F) = 2(-dln7V(> F)/d\nF - 1) is deter- 
mined by the logrithamic slope of the (unlensed) galaxy 
luminosity function N{> F) and is in principle measur- 
able. 

The magnification bias affects both C^j\ through a 

subtle source- lens coupling ijHamanal l2001h . and cfj\ 
However, these impacts are negligible, in the context of 
this paper. The magnification bias has a relatively larger 
effect on c|f'' and modifies Eq. [8] to 



a 



(2). 



-C: 



gG 



gi 



GG 



C: 



gG 



9^C, 



IG 



•gl 1 „ r^GG 



(26) 



tistical error e™'" in shear measurement (A/, 



(a) 



< 



Here gi is the averaged g{F) over galaxies in the i-th 
redshift bin, gj = ja dF)). g[F) is of order unity (e.g. 
IScranton et al.l [20051 ). However, since it changes sign 
from bright end of the luminosity function to the faint 
end. We expect that the averaged gi to be smaller than 
1 for sufficiently deep surveys, gi < 1. 

Our goal is to measure Cf/, with new contaminations 
from the magnification bias. We can apply the same 
weighting of the estimator Eq. [T5| here. On one hand, 
both Cfl and are unchanged by this weighting. On 
the other hand, both C<f ^ and Cf are reduced by virtu- 
ally the same 1 — Q. These behaviors mean that, the es- 
timator Eq. HSj eliminates the combination Cf^ + QiCl^ 
and measures the combination Cf/ -I- giCf^'^ in which the 
term giC^^ contaminates the I-g measurement. 
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giC^'^ can not be eliminated completely, due to vari- 
ous sources of error. An obvious one is the measurement 
error in g{F). At bright end, the galaxy number den- 
sity drops exponentially and lensing modifies N{> F) 
significantly for its steep slope. At faint end, the flux 
measurement no ise is large. Catastrop hic photo-z error 
is also an issue (jSchneider et al.l 120001 ). We will not es- 
timate these errors for realistic surveys. Instead, we ask 
how stringent the requirement on the gi and C^'~^ mea- 
surements should be in order to make the impact of the 
magnification bias negligible. 

Suppose the gi measurement has an error 6gi and the 
Cfl'-^ measurement has an error SC^'~^ , the induced frac- 
tional error in C,^*^ measurement is 
''J 



hi 



< 



= 0(10-3) ( 



^ {\Sg. 



sex. 



Sgi 


+ 




) 


0.1 




10% 





The above relation holds since Cf^'^ < Cf^^ (i < j), 
hi = 0(1) and H/ijA, = 0(10-2). (l) If itself is 
small {\gi\ ^0.1), then there is no need to correct for the 
9i^u'~^ term since its influnece is at the level of 0.1% and 
thus negligble. (2) If gi is large, but it can be measured 
with an accuracy ±0.1, and if O*^'' can be measured with 
10% accuracy, the magnification bias induced error will 
be O(10~3). It can thus be safely neglected, compar- 
ing to the minimum statistical error in Cf^^ (Fig- El) or 
to other residual errors of the self-calibration technique 
(Fig. [3]&|4]). Direct measurement of gi from the observed 
(lensed) fiux galaxy distribution in the redshift bin and 

the approximation C^P ~ O^?'^ likely meet the require- 
ment, unless the II contamination is larger than 10%. (3) 
If the II contamination is larger than 10% and the mea- 
surement of is heavily polluted, a more complicated 
method may work. We can split galaxies into flux bins 
and perform the above analysis to each flux bin. Since 
g{F) changes in a known way across these flux bins, we 
are in principle able to eliminate the giC^'^ term, com- 
bining all these measurements. For example, one can 
flnd an appropriate weighting function W{F) such that 
{g{F)W{F)) = 0. Although this method requires more 
accurate g{F) measurement, it does not require measure- 
ment on Cf^'^ and thus avoids the II contamination and 
other associated errors. 

Based on the above arguments, we expect that our self- 
calibration technique is safe against the magnification 
bias, although extra care is indeed required. 

4.2. Catastrophic photo-z error 

Numerical calculations we perform in this paper only 
consider a Gaussian photo-z PDF. Observationally, the 
photo-z PDF i s more complicated, with non-negligible 
outliers (e.g. iJouvel et al.l l2009t iBernstein fc Hutered 
|2009| ). The existence of this catastrophic error affects 
the self-calibration technique through two ways. (1) It 
affects the accuracy of the Q estimation. (2) It affects 
the scaling relation Eq. 1101 As shown in the appendix 
iRl and further discussed in H'd.'Sl a key condition in de- 
riving Eq. [To] is that the true galaxy distribution in a 



given photo-z bin is sufficiently narrow and smooth. So 
likely catastrophic error leads to degradation of the scal- 
ing relation (Eq. [TU])."' 

From the appendix B, catastrophic error introduces 
bias to Q, mainly through its impact on 77. Stage IV 
lensing projects ne ed to control the o utlier rate /cat to 
~ 0.1% accuracy (jHearin et al.l l2010l ) in order for the 
induced systematical errors to be sub-dominant. If it 
is the case, we are able to perturb the photo-z PDF 
p{z\z^) in Eq. IB51 We choose |z — z^| > A as the 
criteria of the catastrophic error and then have /cat = 

^p{z\z^)dz + J^_^^p{z\z^)dz. Since /cat < 1, 
from Eq. IB5[ we flnd that the induced bias SQ = 0(/cat)- 
As long as the goal |/cat| ^ 0.1% can be achieved, the 
induced error in Q is less than 1% and hence not a signifi- 
cant source of error in the self-calibration. Furthermore, 
we are able to infer the statistically averaged photo-z 
PDF through self- and cross-calibrations of photo-z er- 
rors, even wit h the presence of large catastrophic errors 
(e.g. [Schn eider et all 120061: iNewmanI 120081: iZhang et al] 
120091: iBeniamin et al l |2010f). Since Q can be predicted 
given the photo-z PDF, we are able to reduce the possible 
bias in Q, even if the actual /cat ^ 0.1%. 

The catastrophic error also affects the scaling relation. 
It biases both C^'~^, through the term Wj and Ui in Eq. 
[KU and C^'^, through the term nf in Eq. [^ Part of 
the effect cancels when taking the ratio of the two. The 
residual error is also of the order 0(/cat)- Hence from the 
above order of magnitude estimation, the bias induced by 
catastrophic error is likely sub-dominant to the major 
systematical error Sfj^j^ in the scaling relation. More 
sophisticated analysis is required to robustly quantify its 
impact. 

4.3. Stochastic galaxy bias 

A key assumption in Eq. [TO] is the deterministic galaxy 
bias with respect to the matter distribution. In real- 
ity there is stochasticity in galaxy distributions, which 
can both cause random scatters and systematic shift in 
the scaling relation. Quantifying its impact is beyond 
our capability, since the galaxy stochasticity, the intrin- 
sic alignment and correlation between the two are not 
well understood. For example, the galaxy bias is likely 
correlated with the strength of the intrinsic alignment, 
since both depend on the type of galaxies. Nonethe- 
less, there are hopes to control its impact. (1) The 
galaxy stochasticity c an in principl e be measured (e.g. 
,Pen ^1998; Fan 200ljBoiioli_fcPenl [2008: Zhang 2009) 
and modeled (e.g. iBaldauf et al.|[2009f ). Measurement 
and modeling of the intrinsic align ment can be improved 
too (e.g. H irata & Scljak 2004; lOkumura et~all 120081: 
iSchneider fc Bridla i2009l) . (2^ Recently I Baldauf et all 
(|2009[ ) showed that, by proper weighting and modeling, 
the galaxy stochasticity can be suppressed to 1% level to 
k ~ 1/i/Mpc. Thus there is promise to control the error 
induced by the galaxy stochasticity in the self-calibration 
to be ~ 1% X flj. This error is sub-dominant to other 
systematical errors, especially the one induced by the 
scaling relation inaccuracy ( tj3.3p . 

However, some forms of catastrophic error brin g be tt er m atch 
in the redshift evolution of the integrands of Eq. lAll & IA3I and 
thus can actually improve the accuracy of the scaling relation. 
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4.4. Cosmological uncertainties 

The self-calibration techniques require evaluation of 
Wij in Eq. [10] and Q in Eq. [151 Both evaluations involve 
the cosmology-dependent lensing kernel Wl{zl, zq) oc 
Jlm(l -I- 2;l)xl(1 - Xl/xg)- Fortunately, we do not need 
strong cosmological priors to evaluate i t. flm has al- 
ready been measured to 5% accuracy (jKomatsu et al.l 
I2OIOI ) and will be measured to below 1% accuracy by 
Planck.^ Stage III BAO and supernova surveys will 
meas ure the distance-reds hift relation to 1% accuracy 
(e.g. lAlbrecht et all 120061 ). So if we take these priors, 
uncertainties in cosmology can at most bias the self- 
calibration at percent level accuracy, negligible to the 
identified ~ 10% scaling relation error in §3.31 We need 
further investigation to robustly quantify the impact of 
uncertainties in cosmological parameters. 

5. DISCUSSIONS AND SUMMARY 

We have proposed a self-calibration technique to elimi- 
nate the GI contamination in cosmic shear measurement. 
It contains two original ingredients. (1) This technique is 
able to extract the I-g cross correlation from the galaxy 
density-ellipticity correlation of the same redshift bin in 
the given lensing survey with photo-z measurement. (2) 
It then converts this I-g measurement into a measure 
of the GI correlation through a generic scaling relation. 
The self-calibration technique has only moderate require- 
ment on the photo-z accuracy and results in little loss 
of cosmological information. We have performed sim- 
ple estimation on the performance of this self-calibration 
technique, which suggests that it can either render the 
systematical GI contamination into a negligible statisti- 
cal error, or suppress the GI contamination by a factor 
of ~ 10, whichever is larger. 

The GI self-calibratio n can be combin ed with the 
photo-z self-calibration (jZhang et al.l l2009t ) for a joint 
self-calibration against both the GI contamination and 
the photo-z outliers. This combination does not over-use 
the information in weak lensing surveys. The GI self- 
calibration mainly use the galaxy ellipticity-density cor- 
relation in the same redshift bin. On the other hand, the 



photo-z self-calibration mainly relies on the cross corre- 
lation between galaxy ellipticity-density correlation be- 
tween different redshift bins. 

More robust and self-consistent evaluation of the self- 
calibration (GI and photo-z) performance requires com- 
prehensive analysis of all relevant errors discussed in ^& 
[31 and possibly more, along with realistic model of galaxy 
bias and intrinsic alignment. We expect that our self- 
calibration technique will still work under this more com- 
plicated and more realistic situation, since the method 
to extract the I-g correlation and the scaling relation be- 
tween I-g and I-G are robus t against the complexitie s 
mmentioned abov e. Recentlv. I Joachimi fc Bridld (|2009| ) : 
iKirk et al] (|2010[ ) proposed simultaneous fittings of cos- 
mological parameters, galaxy bias and intrinsic align- 
ment. Our self-calibration technique can be incorporated 
in a similar scheme. 

Our self-calibration technique only uses the shape- 
density and density-density measurements in the same 
redshift bin to calibrate the intrinsic alignment. Lens- 
ing surveys contain more information on the in- 
trinsic alignment, in the shape-shape correlation of 
the same and between different redshift bins, and 
shape-density correlation between different redshift 
bins. The se information h as be ei i inc orporated by 
iJoachimi fc Sc hneider (20 0. 81. I2009D : iJoachimi &: Bridli 
(l2009h: lOkiimura fc Jinei ([2009h : iKirk et al.l ([2010); 
iShi et all (l2oTolv r to calibrate the intrinsic alignment. 
These techniques are complementary to each other and 
shall be combined together for optimal calibration. 
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APPENDIX 
A: THE SCALING RELATION 

We derive the scaling relation (Eq. [TU)) under the Limber approximation. Under this approximation, the 2D GI 
angular cross correlation power spectrum between the i-th and j-th redshift bins is related to the 3D matter-intrinsic 
alignment cross correlation power spectrum A^j(fc, z) by 



Here, 



z I Wjiz)x{z)ni{z)dz 



Wl{zl, ZG)nj{zG)dzG 



(Al) 



(A2) 



As a reminder, ni(z) is the true redshift distribution of galaxies in the i-th redshift bin and Wl{zl, zg) is the lensing 
kernel. The integral limit runs from zero to infinite, to take into account the photo-z errors. On the other hand, the 
2D angular power spectrum between the intrinsic alignment and galaxy number density in the i-th redshift bin is 



^gi 



k = 



nUz)x{z)^dz = b.{i)^^ 



A2 



k = 



-^,z] nj{z)x{z)'^dz . (A3) 
X\z) J dx 



Agj(fc, z) is the 3D galaxy-galaxy intrinsic alignment power spectrum. In the last relation we have adopted a determin- 



istic galaxy bias bg{k,z) with respect to matter distribution and thus Agj{k,z) = bg{k, z)A^j{k, z). bi{i) is defined 
by the above equation. It is the average of bg{k — £/x, z) over the redshift bin. As long as bg{k, z) does not change 
dramatically, we have bi{£) = bg{k = i/xii Zi), to a good approximation. 

In the limit that the ture redshift distribution of galaxies in the i-th redshift bin is narrow, A^^j (^g/) changes slowly 
and can be approximated as A^j(fc = ^|xi^ Zi) {Agj(k = £/xij Zi)). We then have the following approximations, 



" A 9 

-A"* r 



z% WijXi , 



(A4) 



and 



2tt 



Xi_ 

A,; 



(A5) 



The quantity Wij and A^ are already defined by Eq. [TT] & [TU Based on the above two equations, we derive Eq. [TUl 
whose accuracy is quantified in 
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B: EVALUATING THE Q PARAMETER 

To derive the relation between Cf^^\s and Cf^^ , namely, Q = C^^^\s/C^^ , we will begin with the real space angular 
correlation function. We denote the angular correlation function between the shear at Zq and galaxies at as 
w'^^{9; Zq, Zg). Its average over the distribution of galaxies in the i-th redshift bin is 

= T'^^'"^' dz^ /■-.+A^./2 ^^p^c^^^, z^)n^(z^)n^{z^g)dz^dz^ 

Jzi-Azi/2 Jzi-l:^Zil2 
^Zi^lS.Zi jl pZi-\-/S.Zi l2 poo poo 

dzQ / dzg / dzG / dzg[w^^{0;zG,Zg)p{zG\zG)p{zg\Zg)nf{zQ)nf{Zg)] . (Bl) 

lzi-Azi/2 Jzi-Azi/2 JQ JQ 

Here, p{z\z^) is the photo-z PDF. As a reminder, we have normalized such that 

l*Zi — Azi/2 POO 

dz^nf {z^) — I ni{z)dz = 1 . 

i-Azi/2 Jo 

Since 

w'^^i0;zG,Zg) - J{dr,^{e';zL)Sg{e' +0;Zg))WLizL,ZG)dzL , (B2) 

where (• • •) denotes the ensemble average and in practice denotes equivalently the average over 9 (the ergodicity 
assumption), we have 

nZi+Azi/2 pZi-\-Azi/2 pcsD poo poo 

w%%e)^ / dzP / dz^ / dzG / dzg / dzL (B3) 

Jzi~Azi/2 Jzi-Azi/2 Jo Jo Jo 

{Srn{d ;zL)Sg{e +e;Zg))WL{zL,ZG)p{zG\zG)p{zg\Zg)n[{zQ)n[{Zg) 
{6.m{0' ; ZL)6g{d' + 6; Zg))W^{zL)n,{zg 



dzL / dzg 
Jo 

Notice that the lensing kernel Wl{zl, zg) — when zl > zg- The averge over all pairs with Zq < Zg gives the other 
correlation function, 

^?i'\s{0)= J {S„AO';zL)Sg{0' +0;Zg))W,{zLH{zg)dzLdzgr,{zL,Zg) . (B4) 

Here, 

^ ^ _ 2 j^l^^^^ljl dzg Jg'^^^^l^ dz^ dzGWL{zL, ZG)p{zG\zg)p{zg\zP)S{zg, z^)nf{zg)nf [z^) ^^^^ 
!^l~Azi/l dzg !^l^Azl/2 /r dzGWL{zL,ZG)p{zG\zg)p{zg\zP)nf{zP)nf{zP) 



where the selection function S{zQ,Zg) = 1 if < and S{zQ,Zg) — otherwise. The factor 2 comes from the 
relation 

Ii^-aS dzg '^4p(^G\zg)pizg\zP)ni'izP)n[izP) ^ ^ ^^^^ 

Is!-Ai!/2 dzg Ii^^Az!/2 dzPp{zG\zg)p{zg\zP)S{zg,zP)nf {zg)nP {zP) 

The power spectra Cf^^ and Cf^^\s are the Fourier transform of wf^^ and respectively. To evaluate these 

power spectra, we again follow the Limber approximation, which states that the dominant correlation signal comes 
from Zl — Zg. We then have 



tj^ x{zy 



KaAk^ ^,z x{z)W,{z)n,{z)dz , (B7) 



2^ ^ Kig = ^ j W^{z)x{z)n,{z)r]{z, Zg = z)dz . (B8) 

The quantity that we want to evaluate is Q{t) = C'~^^\s{£)/C'~^^{i). Since it is the ratio of the two power spectra, 
in which A^^, Wi and Ui roughly cancel, to the first order, Q ^ rj. The value of rj is determined by the relative 

contribution to C'^^ from pairs with zg < Zg and pairs with zg > Zg . If the two sets have the same contribution, 
7] = 1 and (5 = 1. In the limit ap ^ Az, the contribution from the pair with zg < Zg to C*^^ approaches to that 
of the pair with zg > Zg . So we have 77 — ?• 1 and Q — > 1. In this limiting case, we will be no longer able to use this 
weighting scheme to separate C*^^ and C^^ . On the other hand, if ap ^ Az, the pair with zg < Zg virtually does not 
contribute to C*^^, we will have 77 — > and Q ^ 0, as would happen for spectroscopic redshifts. As long as Q deviates 
significantly from unity, C// can be separated from Cf^^^ . For a LSST-like survey with Az — 0.2 and crp — 0.05(1 + z), 
we numerically evaluate ri{z,Zg) and Q{1)- We find that Q ^ 1/2 (Fig. [5]). The significant deviation of Q from unity 
is the manifestation of relatively large photo-z error ap , across which the lensing efficiency changes dramatically. 
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C: THE STATISTICAL ERROR IN EXTRACTING C/f 



For the convenience, we will work on the pixel space to derive the statistical error in extracting C// from the galaxy 
shape (ellipticity)-density measurement in the i-th photo-z bin. For a given redshift bin, we first pixelize the data into 
sufficiently fine (and uniform) pixels of photo-z and angular position. Each pixel, with label a, has a corresponding 

photo-z and corresponding angular position 9^ ■ Each pixel also has a measured overdensity + and a measured 



"shear" , Ka + la 



Here, the superscript "N" denotes the measurement noise, e.g., the shot noise. In total, there 



are Np pixels. Following the definition of the angular power spectrum, we have 



^ J2 [Sa + S^] + + exp \te- (Oa - Op) 



= 2A^p2 ^ ^ ^N^^ ^AT] J^^-. (^-^ _ ^ g^^ 

Here, Sap = 1 when > Zp and vanishes otherwise. In the limit that Np ^ 1, "^ap ^a/s 
average Sap — 1/2. 
The C^^ measurement error, from Eq. 1151 is 
1 



(CI) 



7V|,/2. Namely, the 



(1-0) 
1 



'^exp 

aP 



ll ■ {da - Op) [{Sa + S^)iKp +Ip+ 4){2Sap - Q) - [1 - Q)SaIp] 



" (1_Q) ^P' J2 K' - [('^"(''^ + ^^p)+Sa i^P +lp+4 )] {"iSap - Q) 



(C2) 



aP 



The last expression has utilized the relation Sap = 1/2 and the fact that the density-intrinsic alignment correlation 
does not depend on the ordering along the line-of-sightof galaxy pairs. The rms error is 



(AC^^')2 = . 



:Ap4 ^ exp 

aPpcr 



exp 



-ie-{ep~9,) {2Sap-Q){2Sp,-Q) 



x{[{5a{np + 4)+S^ {np +Ip + K^ )] [{Spin, + ^^') + {n, + 1. + ^^)]) . 

Here, (■ • •) denotes the ensemble average. To proceed, we adopt a common simplification in the lensing error analysis, 
namely the Wick theorem for 4-point correlation (which holds strictly for Gaussian field) , 

(ABCD) = {AB){CD) + {AC){BD) + (EC) (AD) ; A,B,C,De 6, 6^, k, k^, / . 
Plug the above equation in and keep all non- vanishing terms, we then have 

1 --4 V- 



^ ^' aPpa 



exp 



li ■ i9a - Op) 



exp 



a- 



0^) {2Sap~Q)(2Sp„~Q) 



X [{5a5p){npn,) + -t- {5^5^){{np + Ip){n, + I,)) + {5^ 5'^ ) ))] . (C3) 

Notice that the sum over {SaKp){SpKa-) terms vanishes, resulting from the definition of Q. The sum over {6aHa){5pKp) 
terms also vanishes, due to the mis-match between the Fourier phases (e.g. the one cx Oa + Oa-) and the angular 

dependence of the correlation f unct ions (e.g. {SaKa) — Wgc{Oa — Oa-)) ■ 

Correlation functions in Eq. IC3I depend on the absolute pair separation, but not on the relative pair orientation. 
For example, (SaSp) = {Spda) = Wg{\Oa — Op\), whose Fourier transform is the galaxy angular power spectrum C^^ . 
This allows us to do the sums above analytically and obtain the final expression of the rms error for a single i mode 



1 



3(1-0)2 



^ (jgg^N^GCN 



(C4) 



In the above expression, we have used the fact that noises only correlate at zero-lag {{up ) oc Sp^ and {6^ ) cx Sp^) 
and the following relations. 



^ i'^Sap - Q){2Spa - 0) ^ (1 - 0)' 



a/3pa 



J2(.^Sap - Q){2Spp - g) ^ (1 - Q)2 + i 



aPp 



1 



^ aP 



^i2SaP~Q)^- 



:(i-0)' + i 



In the above relations, we have neglected terms of the order 0{1/Np) and higher, since the number of pixels Np ^ 1. 
Each term in the r.h.s of Eq. IC4I has specific physical meaning and hence deserves brief explanation. 
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• The first term C^ac^^G ^j^g 

cosmic variance arising from the lensing and galaxy density fluctuations. The 
Q dependence drops out, since both C^^^ and C'-^-'|s sample the same cosmic volume and share the identical 
(fractional) cosmic variance from this term. 

(739 (7GG jg familiar term in the cosmic variance of the ordinary galaxy-galaxy lensing power spectrum. However, 
the other familiar term, C^*^'^, does not show up here. This again is caused by the fact that both C*^^^ and 
C'-^-'ls sample the same cosmic volume and the cosmic variances inducing C^'^''^ cancel in the estimator Eq. [151 

• The last term C33-^C^G,N^i + 1/(1 - Q)^] is the contribution from the shot noise in the galaxy distribution and 
random shape shot noise in shear measurement. The Q dependence can be understood as follows. Such error in 
C^^-* has two contributions, SCa from pairs with > Zq and SCb from pairs with < Zq. The total error 
is {SCa + SCb)/'2- Since they come from different pairs, these two errors are uncorrelated {{SCaSC^) — 0), but 
they have the same dispersion {\SCa\'^) = (I^C'sP) = 2C^^'^C"^*^'^. The factor 2 here provides the correct rms 
noise m C(2), which is csf'^C"^'^'^. Clearly the shot noise error in C^^^ \s is SC a- P lug the above relations into 
Eq. [121 we find that the shot noise contribution is indeed the last term in Eq. IC41 Unlike the cosmic variance 
term, which does not rely on Q, the shot noise term blows up when Q ^ 1. This corresponds to the case that 
the photo-z error is too large to provide any useful information and thus we are no longer able to separate the 
Ig contribution form the Gg contribution. 

• The middle term is the cross talk between cosmic variance and shot noise. One can find similar terms in usual 
cross correlation statistical error analysis. 

Interestingly, when (1 - Q)^ = 1/3 and when C^^ < C"^*^, Eq. [C4]reduces to 

{AC^^f~{C3a+2Ca3^^){C^^ + 2C^^'^), w;/ien g - 1 - VX^ = 0.423 . (C5) 
This expression is identical to the usual expression of cross correlation statistical error, expect for the factor 2. 



